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Time-bin encoded photonic qubits are robust against decoherence in birefringent media such 
as optical fibers [T], which has led to many successful quantum communication demonstrations, 
including a Mach-Zehnder interferometer-based system [2], plug-and-play quantum key distribution 
(QKD) OH], differential phase shift (DPS) [5] and coherent one-way (COW) [6] protocols, quantum 
teleportation [7], and elements of quantum repeaters [8]. Despite these successes, time-bin encoding 
is considered impractical for turbulent free-space or multimode fiber channels, because spatial and 
temporal mode distortions hinder the interference required for analysis of time-bin superposition 
states. Here we present an efficient time-bin qubit analyzer, assisted by passive relay optics, that 
matches the spatial mode of the interfering paths. Our system demonstrates interference visibility 
of 89%, even for strongly distorted optical input modes emerging from multimode fiber under 
variable angles of incidence. We measure a high level of entanglement between a time-bin qubit 
and a separate polarization qubit, thereby demonstrating the feasibility of time-bin based quantum 
communication over optical multimodal links. 


Polarization-encoded photons have been preferred 
for most quantum communication experiments in free 
space mis] due to their robustness against atmospheric 
turbulence [14]. However, over long distances or when 
the communicating parties are in motion, it can be chal¬ 
lenging to maintain alignment of polarization reference 
frames. In addition, building polarization preserving op¬ 
tics and telescopes is often difficult and expensive. In 
particular, active or adaptive optics, utilized to compen¬ 
sate angular fluctuations of the link in real time, may 
incur unavoidable polarization fluctuations m, leading 
to practical challenges for polarization-based quantum 
communications M- 

Time-bin encoding is an alternative which has many 
advantages over polarization. Despite this, it has not 
been implemented for quantum communications over a 
long-distance free-space channels. This is mainly due 
to turbulence-induced effects on the wavefront m and 
path m, which make it difficult to perform interfero¬ 
metric measurements necessary for time-bin state pro¬ 
jection onto superposition bases. For instance, Ursin et 
al. |9] measured turbulence-induced angle of incidence 
(AOI) errors of up to 100 //rad, in addition to scintilla¬ 
tion effects. An even larger contribution may come from 
telescope pointing errors between communicating parties, 
which can be as high as ± 1.04 mrad (0.06 °) for a moving 
platform m- Most importantly, these angle variances 
lead to distinguishable paths in unbalanced Michelson 
and Mach-Zehnder interferometers, which form the basis 
for typical time-bin qubit analyzer (TQA) implementa¬ 
tions. 

We analyze the performance reduction of a TQA with 
variable AOI, assuming a Michelson interferometer with 
path difference AL, for an incident Gaussian beam with 
width incident with an angle a. The angle-dependent 


interference visibility follows 

/ AZyQ tan(Q;) 

V(a) = Voe”^ V2cr(l+tan(Q;)) / ^ (2) 

where Vo denotes the system visibility at zero angle (see 
Appendix). For a = 1.49mm and ALq = 0.60m, chosen 
to achieve a clear separation of the time bins given our 
detector timing jitter and channel-induced dispersion, we 
expect the visibility to drop to 0.70 for a = 1.70 mrad and 
Vo = 0.91. Furthermore, the relative phase between the 
two paths is very sensitive to the AOI, with a predicted 
TT-shift per 349 nanorad input angle variation (see Ap¬ 
pendix). The relationship Eq. 0 is verified experimen¬ 
tally with a single-mode beam (see Fig. 1(b)), generated 
by a continuous-wave (CW) laser at 776 nm. As shown 
in Fig. 1(c), the initial interference visibility of = 

0.91 ± 0.01 decreases rapidly with AOI when no correc¬ 
tion optics are implemented, as expected from Eq. Q. 
Next, the same laser beam is sent through a Im-long 
step-index multimode fiber, thereby distorting it into a 
multimodal beam [20| which mimicks the effect of turbu¬ 
lent atmosphere (Eig. 1(d), see OUT] for comparison). 
Despite lengthy and careful alignment we are only able 
to obtain a maximum visibility of = 0.16 ± O.OI, 

which, as shown in Eig. 1(e), drops with AOI. Current 
solutions to this behavior include spatial-mode filtering 
using single-mode optical fibers, which, however, discard 
most of the impinging photons [21]. These observations 
clearly show that, given the expected angular deviations 
reported for free-space quantum channels, it would be 
technically very challenging to achieve a reliable, stable 
and efficient operation of time-bin qubit receiver using 
standard interferometers. 

These interference challenges are overcome by utilizing 



2 


a) 



^ fiber coupler 

U mirror 

1^ 50/50 beam splitter 

[y| piezo actuator 

0 lens 

^ free-space detector 


AOI (degree) 

0 0.05 0.10 0.15 0.20 




0 0.05 0.10 0.15 0.20 

AOI (degree) 


Figure 1. A time-bin qubit analyzer (TQA). a) Schematic diagram of the relay-optics-assisted TQA. Temporal separation 
between paths is set to 2.0 ns (AL = 0.6 m). b) Image of the incident single-mode Gaussian beam, taken by a beam profiling 
camera (WinCamD-UCD12). c) Interference visibility for the single-mode beam with (green circles) and without (black squares) 
relay optics. Measured visibilities are in good agreement with theoretical prediction of Eq. 0 (dashed black line), d) Image 
of the incident multimode beam, generated by a 1 m-long step-index multimode fiber (Thorlabs M43L01). e) Interference 
visibility for the multimode beam with (green circles) and without (black squares) relay optics. Correcting with relay optics 
not only improves performance at higher AOI but is necessary to enable high interference visibility with a multimode beam. 
Uncertainties are smaller than symbol size. 


relay optics in the long arm of the unbalanced Michelson 
interferometer, as depicted in Fig. 1(a). Effectively, the 
relay optics reverse differences in the evolution of the spa¬ 
tial mode over length AL in the long arm (see Appendix). 
With our new TQA design, for a single-mode beam, an 
interference visibility of 0.91 ± 0.01 is obtained, 

which remains constant as the AOI is varied (see Fig. 
1(c)). The improvement is further confirmed by measure¬ 
ments with a multimode beam (Fig. 1(d)) where the high 
visibility of =0.89 ± 0.01 (Fig. 1(e)) demonstrates 
that the interferometer design is robust against highly 
distorted beams. The measured optical throughput of the 
interferometer is 0.74 including multimode fiber-coupling 
efficiency of 0.87, showing that the relay lenses faithfully 
symmetrize the two paths of the interferometer. 

The utility of our interferometer design, as a working 
TQA for multimodal quantum signals, is demonstrated 
with the experimental setup depicted in Fig. 2. A source 
generates polarization-entangled photon pairs centered 
at 776 nm and 842 nm (see Appendix) with initial en¬ 
tanglement visibilities of = 0.979 ± 0.013 and 
= 0.901 ± 0.012. The 776 nm polarization qubit is con¬ 
verted into a time-bin qubit, resulting in the hybrid state 


|^) = T(|h)|E) + |V)|L)), (2) 


where |H) (|V)) and |E) (|L)) denotes the quantum state 
in which the 842 nm photon is in horizontal (vertical) 
polarization mode and the 776 nm photon in early (late) 
temporal mode (see Appendix). 

The 776 nm photon is sent through a 1 m-long step- 
index multimode fiber to artificially distort the spatial 
mode (see Fig. 3(a)) and the temporal mode, with a 
measured dispersion of about 50 ps, drastically exceeding 
the photon’s coherence time of 3.2 ps [22], prior to enter¬ 
ing our multimode time-bin qubit analyzer (MM-TQA). 
The output of the interferometer is coupled into a mul¬ 
timode fiber for delivery of the photons to the detector 
(Si-avalanche photodiode). 

Verification of the MM-TQA performance is done by 
measuring entanglement visibilities between the polar¬ 
ization and time-bin qubits. Each qubit is first projected 
onto basis states, i.e. |cr±z)(c^±z|: where |cr+z(_ 2 )) denotes 
|H(V)) and |E(L)) for a 842 and a 776nm photon, re¬ 
spectively. The coincidence counts are used to calculate 
a visibilty, defined as V±z = (N±z±z - Nzpz±z)/(N±z±z + 
Nzpz±z), and the average value of Vz = (V+z + V_z)/2 = 
0.952 ± 0.011 is observed. Here, Ny denotes the joint- 
detection counts between a 842 nm photon projected onto 
|cri)(cri| and a 776 nm photon onto |crj)(crj|, where i, j G 
{+z, -z} (see Fig. 3(b)). 

More important for the MM-TQA operation is the pro¬ 
jection of the time-bin qubits onto superposition states. 
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Figure 2. Experimental setup. The polarization-entangled photon-pair source (EPS) is described in the Appendix. For 
projection measurements, a 842 nm photon is directed to a polarization-qubit analyzer (PQA), consisting of a quarter-wave 
plate (QWP), a half-wave plate (HWP), a polarizing beam splitter (PBS), and silicon avalanche photodiodes (Si-APDs) with 
timing jitter of 600 ps. After reflection at a dichroic mirror (DM), by a flip mirror (FM), a 776nm photon is sent either to a 
PQA or a time-bin qubit converter (TQC) followed by a multimode time-bin qubit analyzer (MM-TQA). The purpose of the 
PQA is to measure reference entanglement visibilities with the source of polarization entanglement, which will be compared 
to visibilities with polarization-time entanglement. The TQC maps the 776 nm polarization qubit onto a time-bin qubit (see 
Appendix). For projection measurements, the time-bin qubit, spatially and temporally distorted by a multimode fiber (MMF), 
is then sent to a MM-TQA whose output is coupled into a MMF before detection by a Si-APD. A polarizer (POL) removes 
any possible path-dependent polarization distinguishability. All detection signals are sent to a time tagger for data analysis. 


i.e. \cj±^){cj±^\, where \ct±^) = ^(k+z) ± e*'^|cr_z)). 
To measure the visibility, we vary the relative phase (p' 
between basis states using motorized wave-plates act¬ 
ing on the polarization qubit. A complete scan of the 
phase along the xy-plane of the Bloch sphere is per¬ 
formed within 10 seconds, yielding the average visibility 
Vxy = (V+xy + V-xy)/2 = 0.804 ± 0.006 (see Fig. 3(c)). 

The entanglement verification measurements are also 
carried out for different AOI, demonstrating that the 
measured visibilities (in Fig. 3(d)) are constant within 
experimental errors. Due to the coupling geometry of the 
photons into the multimode detector fiber, the photon- 
collection efficiency of 0.87 decreases with AOI and drops 
to zero at ± 0.24 degrees. However, the fact that the in¬ 
terference visibility is essentially invariant to variation of 
AOI confirms the robustness of our MM-TQA. 

Without relay optics, a variance of AOI also introduces 
phase fluctuations of the time-bin qubit analyzer. For 
instance, Dixon et al. [23] measured an angular deflec¬ 
tion of 400 femtorad using interferometric phase mea¬ 


surements. From our theoretical model, we anticipate 
5 7 r-shift with an AOI of only 1.75 /irad ( 0 . 0001 °) (see 
Appendix). In order to assess the phase stability of our 
MM-TQA with AOI, we measure expectation values, de¬ 
fined as = (N+ 0 / 0 -N_ 0 / 0 )/(N+ 0 / 0 +N_ 0 / 0 ), for AOIs 
changing from - 0 . 20 ° to + 0 . 20 ° continuously over 20 sec¬ 
onds. As shown in Fig. 4, the measured values re¬ 
main almost constant within experimental errors, show¬ 
ing that our MM-TQA prevents the AOI-caused phase 
fluctuations. Moreover, by performing continuous pro¬ 
jection measurements over a half an hour, we also test the 
long-term stability of our passively stabilized MM-TQA 
and observed the measured F^-values change very slowly 
during the measurement period (see Supplementary In¬ 
formation (SI)). Hence, this MM-TQA could be used for 
reference-frame-independent (RFI) QKD in which secret 
keys are extracted under slow drift of measurement bases 
121H27|. 

The suitability of the multimode TQA for quantum 
communications is further substantiated by examining 
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Figure 3. Experimental results for entangled photons analyzed with the novel MM-TQA. a) Spatial mode of 776 nm 
photons before entering the MM-TQA. The image is captured by a Hamamatsu electron multiplier (EM) CCD camera (C9100- 
13). b) Joint detections for the projection a±z <^cr±z as a function of detection-time difference between a 776 nm and a 842 nm 
photon, c) Joint detections for the projection (g) cr^^ as a function of phase 0' of a polarization qubit. Visibilities V±xy are 
obtained from sinusoidal fittings. Single counts remain essentially constant as we scan the phase, d) The measurement b) and 
c) are repeated for different AOIs. Blue diamonds are obtained from projection measurements along the computational basis 
(z-axis of Bloch sphere, as seen in (b)) and red circles superposition basis (xy-plane of Bloch sphere, as seen in (c)). Dashed 
blue and red lines are reference visibilities measured directly with the source of polarization-polarization entanglement. We 
maintain high entanglement visibility (close to source visibility) despite the high multimode nature of the incoming photons. 


a CHSH-Bell inequality [28] from measurements on the 
hybrid entangled state. We prepare two bases on xz- 
plane of a Bloch sphere for the polarization qubit, i.e., 
|<^z±x)(<^z±x|- For the time-bin qubit, one basis {\(Jz){o'z\) 
is set deterministically, while the second one (|crx)(crx|) is 
made to drift using a heat gun, owing to the absence of 
active phase control in our current TQA. A search for the 
maximally achievable CHSH-Bell parameter is performed 
by scanning the second basis {\(7cj)){( t( f)\) along the xy- 
plane of the time-bin states. A maximum value of iSexp 
= 2.42 ± 0.05 is found. This is in good agreement with 
the predicted parameter 5'theo = 2.47 ± 0.02 calculated 
assuming the state described in Eq. is exposed to 
noise modeled by an asymmetric depolarization channel 
(see SI). 

A more rigorous verification of entanglement takes 
practical assumptions into account, including the photon 
loss in the channel and path-dependent transmission in 
the MM-TQA. We numerically find minimum visibilities 
required to detect the presence of entanglement for arbi¬ 
trary 2 X 3-dimensional quantum system (see SI). Mea¬ 
sured visibilities Vz = 0.95 ± 0.02 and Vxy = 0.81 ± 0.02 
averaged over various AOIs are well above the obtained 
classical bound, revealing clear evidence of entanglement 
(see Fig. 5). Note that the obtained numerical results are 
valid regardless of the system efficiency and mismatched 
transmission between the two paths of TQA. 

While theoretically Vo ~ 1, the performance of the 
current MM-TQA shows about Vo = 0.9, which is ex¬ 
pected to be improved considering several approaches. 
First, the overlap of the spatial modes in the interfer¬ 


ometer can be improved with careful custom design and 
selection of optical elements and optimization of beam 
diameters. Second, appropriate mode matching optics in 
both paths will improve the interference quality because 
of dispersion symmetrization at the expense of increasing 
complexity of the system. 

In conclusion, we demonstrated a novel time-bin qubit 
analyzer, robust against multimode distortions on the 
quantum link as expected from turbulence-induced ef¬ 
fects, i.e. wavefront distortions and path-delay fluctu¬ 
ations, as well as telescope pointing errors. The main 
benefit of our analyzer is that it is entirely based on 
passive optics without any active feedback systems for 
beam correction. Despite strong spatial and temporal 
distortions of the incident photon modes, we observe a 
robust entanglement visibility of 0.80 in the superposi¬ 
tion bases which remains constant with variable angles 
of incidence. At the same time, our MM-TQA has a 
very high throughput efficiency of 0.74 from input to out¬ 
put, mainly limited by optical surface losses which can 
be improved. Our results open the door for implement¬ 
ing time-bin based quantum communication experiments 
over turbulent free-space and multimodal channels such 
as with moving platforms, including aircraft and satel¬ 
lites, or through non-polarization maintaining windows. 
Furthermore, as a quantum state receiver, the multimode 
time-bin qubit analyzer should also be suitable for COW, 
DPS, and reference-frame-independent QKD protocols. 
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Figure 4. Phase stability of MM-TQA Expectation values 
(green circles) are measured as the AOI is continuously varied 
from -0.2° to +0.2° over 20 seconds. The inset shows the ex¬ 
pectation value without relay optics as a function of AOI. Due 
to AOI-induced phase fluctuations, the value rapidly changes 
with AOI and thus yields an average value of zero (red circle). 
Measured visibility Vxy = 0.80 bounds the value Ecfy. These 
phase fluctuations are corrected by the relay optics, allowing 
a near constant expectation value. 


APPENDIX 


Analysis of an uncorrected time-bin qubit analyzer. 

Let us consider an unbalanced Michelson interferometer 
with long path Li and short path I/2. The path differ¬ 
ence for normal incidence is ALq = 2(I/i — L 2 ). Simple 
ray tracing shows for an incident angle a 7^ 0 the paths 
through the interferometer change with a according to 


AL{a) = 


ALn 


+ 


1 — tan(a) 


cos((a) cos((a) + sin((a) 
+ 5{a) tan ~ ^ 


where 


6{a) 


ALq tan((a) 
1 + tan((a) 


is the lateral offset between the two rays coming from 
the short and long arms at the output port of the beam 
splitter. The angle-dependent path difference affects the 
photon counts in the phase-dependent basis, as the inset 
in Fig. 4 shows. The lateral offset reduces the interfer¬ 
ence visibility, as shown in Fig. 1. We can see this if the 
incoming ray is replaced by a single-mode Gaussian beam 
with an amplitude of a and variance The intensity 
at the output is then given by 




2 2 
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Figure 5. Entanglement verification. Negative Partial 
Transpose (NPT) criterion [29l [30] is used to obtain the re¬ 
quired entanglement visibilities, certifying the presence of en¬ 
tanglement in arbitrary 2 x 3-dimensional quantum state. 
Our experimental results for various angles (green circles) are 
well above the classical bound (black line). 


from which we obtain the visibility in Eq. Q. Here, (f) 
denotes the relative phase between the paths. 

Principle of relay optics. The relay-lens system, con¬ 
sisting of free space (FS) and lens (L) transmission. 





FS 




FS 


= IL 


can be inserted into the long arm of an unbalanced 
Michelson interferometer (see Fig. 1(a)) in order to 
remove the spatial evolution caused by the path-length 
difference. Here, / denotes the focal length of the lenses. 
Note that image inversion at mirrors is excluded as 
this happens in both paths with no effect on the path 
indistinguishability. 


Source of polarization-entangled photons. Light from 
a 404 nm CW laser with an average power of 6mW 
pumps a periodically poled potassium titanyl phosphate 
(PPKTP) crystal inside a Sagnac interferometer. This 
generates polarization entangled photon pairs at 776 nm 
and 842nm via type-H SPDC: = ^(|HV) + |VH)), 
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where H(V) denotes horizontal {vertical) polarization. 
The downconverted photons are spectrally filtered with 
12 nm band-pass filters (BPF) and pump photons are 
removed with long-pass filters (LPFs). The ratio of 
coincidence-to-single counts is measured to be 0.07. 

Time-bin qubit conversion. The 776 nm photon from 
the polarization-entangled photon pair is sent through a 
single-mode fiber (SMF) to the time-bin qubit converter 
(TQC). A polarization controller (FPC) ensures the 
faithful transmission of the polarization state. At the 
PBS of the TQC, the photon is either transmitted or 
reflected into the short or long path, respectively. To 
ensure both photons leave the interferometer at the 
desired output port, a quarter-wave plate is inserted 
in each path. Directly after the interferometer the 
photon passes through a polarizer (POL) set to an equal 
superposition between the polarizations. This allows, at 
the cost of 50% transmission loss, to erase polarization 
information for each time-bin state, which completes 
the map |H) i-^ |L) and |V) |E) from polarization to 

time-bin state. The time-bin qubit is then coupled into 
a SMF with a collection efficiency of 0.67, which gives a 
total transmission through the TQC of 0.24. 
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SUPPLEMENTARY INFORMATION 

A. Estimation of CHSH-Bell parameter 

The non-classicality of time-polarization entanglement 
is bounded with an estimate of the Bell-CHSH inequal¬ 
ity violation [T]. Despite the absence of an active phase 
control, required to set measurement bases for the time- 
bin qubit deterministically, we search for the maximally 
obtainable violation by varying the measurement basis. 
More specifically, we first set the measurement basis for 
the polarization qubit to Ai = |crz+x)(<^z+x| using wave 
plates, and slowly and continuously change the path- 
length difference of the multimode time-bin qubit ana¬ 
lyzer (MM-TQA) by externally heating it. This allows us 
to scan projection measurements for the time-bin qubit 
in the superposition bases. Fig. 6 (a) shows coincidences 
between a polarization qubit (two detectors, i.e. D1 and 
D 2 ) and a time-bin qubit (three temporal modes). De¬ 
tections in the early/late bin {middle bin) corresponds to 
a projection of the time-bin qubit onto Bi = \(j/i{(j/\ {B 2 
= |( 70 )(cr 0 |). Owing to the absence of the second output 
of the TQA, we consider all possible expectation values 
E{Ai,B/ between any two points in time, i.e. ti and ^2 
(see Fig. 6 (b)), which is defined as 


N++ + N:t- - N+- - N:t+ 

-T- 

^ N+++ -h N+“-1-N“+ 


( 3 ) 


Here, Ny are the coincidence counts for the projections 
Ai (g) Hj, where i,j g{1,2} and superscript(-h, —) denotes 
two outcomes of the projection measurement. Among 
all the computed expectation values, we find the abso¬ 
lute maximum expectation value. We then change the 
measurement basis for the polarization qubit to A 2 = 
|crz-x)(<^z-x| and repeat the procedure. Finally, we com¬ 
pute the CHSH-Bell inequality parameter 


5 = |^(Ai, Hi) - H(Ai, H 2 ) + H(A2, Hi) + H(A2, H 2 )| 

( 4 ) 

and find the value of Hexp = 2.42 ± 0.05, which is clearly 
above classical bound S = 2. To see whether this value 
agrees with the measured visibilities, we model the two- 
qubit state with noise, described by an asymmetric de¬ 
polarization channel, on a time-bin qubit. The output 
state is then described by 


Pout — (1- Pimii’i+ Yi Pj(10(7j)|?/))(tA|(10crj), 
j=x,y,z j=x,y,z 

( 5 ) 


where pj(j=x,y,z) and ctj are the depolarization proba¬ 
bility and single-qubit Pauli operator along the channel 
j, respectively. Here, \/) denotes the input state, de¬ 
scribed in Eq. (§ of the main text. Assuming unbiased 
depolarizations in the superposition bases, i.e. (px = 
Py) = pxy, we calculate expectation values H(Ai,Hj) = 
Tr (Pout Ai(g)Hj) for given measurement bases. Using the 
definition of visibility, we further represent the CHSH- 
Bell parameter Htheo = V^(Vz + Vxy) as a function of 
entanglement visibilities. We find Htheo = 2.47 ± 0.02 is 
in accordance with our measured value Hexp = 2.42±0.05. 


B. Long-term phase stability of the MM-TQA 

Our multimode time-bin qubit analyzer is passively 
stabilized by enclosing it with black hardboard. In order 
to assess the phase stability of the MM-TQA, we per¬ 
form joint-projection measurement onto the superposi¬ 
tion bases over a half an hour. The time-bin qubit is pro¬ 
jected onto a \(J(f)){(J(f)\ and the polarization qubit alter¬ 
natively between or |cr 0 /+^/ 2 )(cr 0 /+ 7 ,/ 2 |. This 

allows US to calculate the expectation values, defined as 
^(f)',(f) — {'^(f) —(j)'(f)^/{'^-\-<p'(j) T where Ny 

denotes joint-detection counts when a 842 nm photon is 
projected onto |cri)(( 7 i| and a 776 nm photon onto |. 

As shown in Fig. 7, the average expectation value 
\j ^ 0 '+ 7 r /2 0 always higher than 0.65, 

which is well above the required value Vxy for verifying 
entanglement, given entanglement visibility Vz = 0.952 
(see Fig. 5 in the main text). 


C. Entanglement verification 

As entanglement is fundamentally important and an 
effective distribution of entanglement is a necessary con¬ 
dition for secure QKD [ 2 ], we need to verify whether or 
not entanglement is present in our experiment. This is 
especially important in the absence of a complete security 
analysis of a QKD implementation. We assume that the 
spontaneous parametric down-conversion process gener¬ 
ates a pair of photons with negligible multiple-photon- 
pairs event. Each photon is a polarization qubit and the 
pair of photons is potentially entangled. By detecting one 
photon in the pair, we can herald the generation of the 
second photon. After the conversion from polarization 
qubit to time-bin qubit, the second one is transmitted to 
the MM-TQA. Suppose that Alice holds the polarization 
qubit while Bob holds the time-bin qubit. To include 
conversion and transmission losses of the time-bin qubit, 
we enlarge the dimension of Bob’s system from 2 to 3 by 
adding a dimension corresponding to no photon arriving 
at Bob. Hence, the final state p shared by Alice and Bob 
is a 2 X 3-dimensional state. We need to verify whether or 
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Figure 6. Estimation of the CHSH-Bell parameter, (a) Yellow and light blue lines are coincidences from joint projections 
<7z+x 0 a^z{early temporal bin), green and orange lines from 0 a-z{late temporal bin), and purple and blue lines from 
(7z+x <S> cFcf){middle temporal bin). Black dashed lines are times at which maximal expectation values are extracted, (b) Surface 
plot of calculated expectation values for the projections in (a), c) Yellow and light blue lines are coincidences from joint 
projection az-x 0 cr+z, green and orange lines from ctz-x 0 cr_z, and purple and blue lines from ctz-x (8) Black dashed lines are 
times at which maximal expectation values are extracted, (d) Surface plot of calculated expectation values for the projection 
measurements in (c). The Beh-CHSH parameter is calculated using the maximum expectation values. The measurement 
duration is chosen arbitrarily and yields a violation of the inequality of 2.42 ± 0.05. 


not the state p is entangled using only the measurement 
results Vz = 0.952 ± 0.011 and Vxy = 0.804 ± 0.006 with¬ 
out further assumptions on the state. Since the measure¬ 
ments of Alice and Bob are block-diagonal with respect 
to the subspaces of total photon number, as we will show 
below in Eq. and Q, we can also assume without loss 
of generality that the state p shows the same structure. 
This follows from the fact that the measurement struc¬ 
ture allows us to assume that a quantum non-demolition 


measurement of the total photon number is executed be¬ 
fore the actual measurement itself. 

In order to verify entanglement, we need to know how 
to accurately describe the measurements on the polariza¬ 
tion qubit and on the time-bin qubit. For the polarization 
qubit, we measure it in the horizontal/vertical or diag¬ 
onal/anti-diagonal basis, i.e. along the z- or x-axis in 
the Bloch sphere. In the horizontal/vertical basis, these 
measurements are represented as 


"1 O' 


"0 o' 

1 

'1 1' 

1 

■ 1 -1' 

0 0 

,Mv = 

0 1 

,Mo = - 

1 1 

, and Ma = - 

-1 1 


( 6 ) 
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Figure 7. Long-term phase stability of our multimode time-bin qubit analyzer. Red squares are expectation values 
for projection measurements ^< 74 ) and blue hexagons for cr 0 /+ 7 r /2 0 Each measurement is averaged 

over 3 minutes. Black circles are average expectation values The phase drifts slowly, on the order of 7 r /2 

over an half an hour, showing the stability of our MM-TQA. The average expectation value is always higher than the required 
expectation value for entanglement verification. 


where the subscript indicates measurement outcome, and 
H, V, D, or A denotes the horizontal^ vertical^ diagonal^ 
or anti-diagonal polarization. On the other side, for the 
time-bin qubit, the photon loss in the long path or the 
short path of the MM-TQA could be different from each 
other. Hence, the operators corresponding to measure¬ 
ment of the time-bin qubit in the early / late-i\m.e basis or 


in the superposition bases could deviate from the ideal 
case. Without loss of generality, we can choose the rel¬ 
ative phase between the early-t\m.e and late-t\m.e basis 
states in the superposition basis to be zero. Therefore, 
in the basis where the basis states are no photon, one 
photon in the early bin and one photon in the late bin, 
these measurements can be written as 


Me 


M0 


"0 

0 

0 “ 

II 

"0 

0 

o' 

0 

hs 

0 

0 

0 

0 

_0 

0 

0 _ 


_0 

0 

Vi_ 


^3x3 — Me — Mi^ — Mx, 


where the subscript E, L, X, or 0 means that the mea¬ 
surement outcome is early time, late time, the superpo¬ 
sition of the early time and late time, or no detection, 
respectively. No-detection events are due to detection 
inefficiency and the absence of the second output in the 
MM-TQA. In Eq. 0, rji or rfs is the respective trans¬ 
mission efficiency in the long path or the short path of 
the MM-TQA. Note that in our experiment rji and rjg are 
very close to each other. 



and 


(7) 


After knowing the description of Alice’s and Bob’s joint 
state p and also that of their measurements, we can verify 
entanglement by the negative partial-transpose (NPT) 
criterion [3]. The NPT criterion is used because this 
criterion is satisfied if and only if a state is entangled, 
given the state is 2 x 2- or 2 x 3-dimensional [4]. The 
NPT criterion has been applied to verify entanglement 
in QKD systems, such as in [5] . Explicitly, we verify 
entanglement by solving the following semi-definite pro¬ 
gram (SDP): 
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find 

subject to 


P 

p>0 

Trip) = 1 
> 0 

Tr[p(MH (8) M'e — My (8) M-^)] = V+z Tr[p(M'H (8) Me + My (8) Me)] 
Tr[p(Mv (8) Ml - Mh (8) Ml)] = V_z Tr[p(Mv (8) Ml + Mr (8) Ml)] 
Tr[p(MD (8) Mx — Ma (8) Mx)] = Vxy Tr[p(MD (8) Mx + Ma (8) Mx)], 


( 8 ) 


where F is the partial-transpose operation on a subsys¬ 
tem, such as on the polarization-qubit subsystem, and (8) 
denotes the tensor product. Note that, we formulate the 
last three constraints according to the measured visibil¬ 
ities. The first two of them are based on entanglement 
visibilities V+z and V_z conditioned on measurement out¬ 
comes of the time-bin qubit being early time and late 
time, respectively. The last constraint is based on entan¬ 
glement visibility Vxy, where the time-bin qubit comes 
out in the middle bin. Since the MM-TQA has only one 
output, we cannot differentiate the case when the photon 
comes out from the second output if this output exists 
from the case when the photon is lost over the transmis¬ 
sion. Hence, we cannot formulate two constraints based 
on Vxy. 

In our experiment, we verified that within experimen¬ 
tal errors the visibilities V+z = V_z. So, for solving the 
SDP program in Eq. we set V+z = V_z = Vz- Using 
the measured results Vz = 0.952 and Vxy = 0.804, the 
SDP program in Eq. (§ is not feasible, signifying that 
the state p must be entangled. Eurthermore, by numeri¬ 
cally checking over which values of Vz and Vxy the SDP 
program in Eq. is not feasible, we are able to up¬ 
per bound the required visibilities Vz and Vxy in order to 
certify the presence of entanglement in the system. The 
numerical results are shown in Eig. 5 of the main text. 
Erom this figure, one can see that our visibility result 
at any observed incident angle witnesses entanglement 
with high confidence. Einally, we would like to note two 


points. Eirst, the constrains considered in Eq. (§ are 
independent of the transmission or conversion loss of the 
photon arriving at the MM-TQA, and even independent 
of the common photon loss in the two different paths of 
the MM-TQA. Therefore, the upper bounds on the vis¬ 
ibilities Vz and Vxy obtained for verifying entanglement 
are independent of all of these different losses. Second, 
we numerically verified that the classical boundary, as 
shown in Eig. 5 of the main text, is even independent of 
the relative loss between the two paths of the MM-TQA. 
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